An Ultracold Strongly Coupled Gas: a Near-ideal Liquid 
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Feshbach resonances of trapped ultracold alkali atoms allow to vary the atomic scattering length 
a. At very large values of a the system enters an universal strongly coupled regime in which its 
properties — the ground state energy, pressure etc. — become independent of a. We discuss transport 
properties of such systems. In particular, the universality arguments imply that the shear viscosity of 
ultracold Fermi atoms at the Feschbach resonance is proportional to the particle number density n, 
and the Plank constant %: rj = hna v , where a v is a universal constant. Using Heisenberg uncertainty 
principle and Einstein's relation between diffusion and viscosity we argue that the viscosity has the 
lower bound given by a n < (67T)" 1 . We relate the damping of low- frequency density oscillations of 
ultracold optically trapped 6 Li atoms to viscosity and find that the value of the coefficient a v is 
about 0.3. We also show that such a small viscosity can not be explained by kinetic theory based 
on binary scattering. We conclude that the system of ultracold atoms near the Feshbach resonance 
is a near-ideal liquid. 



I. INTRODUCTION 

A presence of the Feshbach resonance in binary colli- 
sions of alkali atoms at ultracold temperatures provides 
an experimental control of an effective strength of the 
inter-particle interaction. In the case of ultracold dilute 
atomic systems the effective interaction strength is given 
by an s-wave scattering length a. For a gas a "diluteness 
parameter" na 3 , where n is the number density, is much 
less then unity.* An atomic system can be tuned to a 
regime of the Feshbach resonance using an external mag- 
netic field. On resonance the scattering length and "di- 
luteness parameter" formally diverge. The system can no 
longer be described as a dilute weakly interacting Fermi 
or Bose gas. This regime is often referred to as strongly 
coupled regime. 

As external magnetic field is tuned to a critical value 
corresponding to the Feshbach resonance the effective di- 
mensionless coupling constant describing two-body inter- 
action formally diverges. In the case of ultracold alkali 
atoms this coupling constant g is given in terms of the s- 
wave scattering length a by g = 4irh 2 a/m, where m is the 
mass of an atom. Since the coupling constant g goes to 
infinity at the Feshbach resonance it is natural to expect 
that the thermodynamic properties of the atomic system 
become independent of g and consequently of the nature 
of the two-body interaction [1]. Instead, these proper- 
ties are given in terms of the universal constants times 
the appropriate dimensional factors which for ultracold 
atoms are given in terms of the remaining parameters, 
namely, the Plank constant fi, particle density n, and an 
atomic mass m. 

For example, the pressure p(n) of ultracold atoms at 
the Feshbach resonance is proportional to that of the 
ideal Fermi gas, Po(n) ~ n 5 ^ 3 /m, 
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*A typical density of ultracold trapped atomic gases is 
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where j3 is an universal dimcnsionless constant. For ultra- 
cold fermionic atoms (3 « 0.5. The "universality" princi- 
ple means that the value of (3 is the same for any strongly 
coupled system of fermions, such as cold atoms or dilute 
neutron matter. This constant is expected to depend on 
the particle spin. 

An atomic system near the Feshbach resonance is an 
example of a quantum many-body system with a large or 
formally infinite dimensionless parameter describing two- 
body interaction. Other examples include: a strongly 
coupled quark-gluon plasma (sQGP) [5], a strongly cou- 
pled classical electromagnetic plasma and a strongly cou- 
pled Af = 4 supersymmetric Yang-Mills gauge theory. 
These systems are briefly discussed in section V. 

In all of the above systems one finds that the dif- 
ference between weak and strong coupling regimes in 
the total (free) energy, although nonzero, is neverthe- 
less much less dramatic than the corresponding differ- 
ence in their transport parameters — viscosity, diffusion, 
etc. — which change by orders of magnitude. In these 
cases the strongly coupled matter behaves as a liquid, 
with the the mean free path comparable or even smaller 
than the inter-particle distance. 

The goal of this paper is to show that the transport 
properties of a system of ultracold atoms change dra- 
matically as the system is driven from a weakly inter- 
acting regime to the regime of the Feshbach resonance. 
As we will discuss in details this change is manifested 
in a steep decrease by about two orders of magnitude of 
the damping rate of the collective oscillations of a cloud 
of 6 Li atoms near the Feshbach resonance observed in 
the experiments of Kinast et al. and Bartenstein et al. 
This decrease cannot be attributed to the phase tran- 
sition that the system is expected to undergo at a low 
enough temperature. The maximum of the oscillations 
is near the surface of the atomic cloud where the tem- 
perature is above the local phase transition temperature. 
Furthermore, the transport coefficients such as viscosity 
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and diffusion display universal properties near the Fesh- 
bach similar to the pressure and energy density. 



II. THE SHEAR VISCOSITY OF ULTRACOLD 
FERMIONIC ATOMS FROM OSCILLATION 
DAMPING 

In this section we use the equations of hydrodynam- 
ics to describe the damping rate of the collective exci- 
tations of the ultracold Fermi atoms near the Feshbach 
resonance. A spectacular manifestation of the hydrody- 
namical behavior of ultracold optically trapped atomic 
clouds near the Feshbach resonance is the elliptic flow 
observed after the trap is switched off [6,7]. Such flow is 
also observed at higher temperatures. Thus it is not by it- 
self related to quantum low-temperature phenomena like 
condensation or superfluidity but rather it is a manifes- 
tation of a strongly coupled regime. A good quantitative 
description of the elliptic flow was obtained in a frame- 
work of ideal hydrodynamics (i.e. neglecting viscosity). 
In principle, this comparison by itself may provide an 
upper limit on the viscosity. 

Additional evidence for the hydrodynamical behavior 
is provided by the experiments of Kinast et al. [2,3] and 
Bartenstein et al. [4] on collective oscillation of ultracold 
trapped 6 Li atoms in the vicinity of the Feshbach res- 
onance. The measured frequencies of the lowest modes 
are in a good agreement with the predictions of ideal 
non-suprefluid hydrodynamics. 

Our main point is that the damping of these collective 
excitations can be described by including the dissipative 
(in particular, viscous) terms in the hydrodynamic equa- 
tions. As will be shown in the next section, the viscous 
terms are very small at least for the low frequency col- 
lective excitations. Thus, these terms can be included 
perturbatively. 

If the collision rate of atoms is large enough to establish 
a local equilibrium the collective vibrations of the atomic 
cloud can be described using standard hydrodynamical 
theory [9] . The collective vibrations are described by the 
local density n(r), pressure p(r) and velocity v(f) which 
are the solutions of the continuity equation, Euler equa- 
tions of motions and the equation of state, respectively, 

dn 



m— +V - (ranv) = 0, 
ran — + ran (v ■ V) v - 



P 



= An' , + 1 . 



(2) 



where V = (1/2) to ^ ojf rf is the harmonic trapping 
potential, A is a constant, and 7 is the polytropic in- 
dex. As discussed in [1] a cold gas of strongly interacting 
fermions at the has a universal equation with polytropic 
index 7 = 2/3. 

The lowest collective modes correspond to the small 
vibrations of the density, pressure and velocity around 



their equilibrium values, n eqi p eq and v eq = 0. These 
values are determined by the static limit of the Eqs. (2) in 
which the Euler equation takes the form Wp eq + n eq VV = 
0. The equilibrium density is 
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and zero outside. In equation (3) n eq (0) is the equilib- 
rium density of the atomic cloud at the center of the 
harmonic trap (r = 0) and 
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are the radii of the cloud with p eq (0) being the equilib- 
rium central value of the pressure. These radii can be 
expressed in terms of the chemical potential fi. Using 
the Gibbs-Duham relation, dp = ndfi (valid at constant 
temperature) and the equation of state, p = An 1+1 , one 
can easily show that 
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so that the Thomas-Fermi radii are 
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independent of 7. 

The chemical potential ji in Eqs. (6) and (7) includes 
interactions of the trapped atoms in the vicinity of the 
Feshbach resonance. The effect of the interaction in the 
unitary limit is to scale the chemical potential of the non- 
interacting trapped Fermi gas [6,8]: 



(8) 



where (3 w 0.5 is the universal constant discussed in the 
introduction and //°) is the chemical potential of TV spin- 
1/2 noninteracting fermions confined in a harmonic trap 
[10], 



u (Q) = hu> (3-ZV) 1 ^ 3 = ks Tj 



F , 



(9) 



where tD 3 = cji(jj 2 ^3 and ks is the Boltzmann constant. 
The last equation also defines the Fermi temperature Tp 
of the system. 

Linearizing Eqs. (2) to describe small density n = 
n eq + Sn e wt and velocity v e wt oscillations one can ob- 
tain the following equations for the density and velocity 
amplitudes: 
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The form of the collective vibrational modes which are 
the solutions of Eq. (10) depend on the symmetry of 
the trap potential. If the confining potential is isotropic, 
lui = lu 2 — ^3, the collective modes have spherical sym- 
metry and can be characterized by definite angular mo- 
mentum and its z component, I and m. The monopole 
mode, I = m = 0, has a velocity profile which is pro- 
portional to r. Such modes are referred to as breathing 
modes. The dipole mode I = 1 involves the motion of the 
center of mass of the cloud and is not usually excited in 
the experiments we are interested in. There are five de- 
generate quadrupole modes corresponding to I = 2 with 
m = 0,±1,±2. 

We will discuss the collective modes that are excited 
in an axially symmetric trap with u>i = lu 2 = w r and 
u>3 = u) z = Aov, where A is a constant. In the experi- 
ment of Bartenstein et al. [4] A is 0.03^. In such traps 
the collective modes corresponding to different angular 
momenta but the same z-component m are mixed. 

To find these modes we look for the solutions of 
Eq. (10) in the form* 



v = (a x x, a y y, a z z) = (a r x, a r y, a z z) . 



(11) 



The set of equations in Eq. (10) for v reduces to the 
secular equation for the eigenfrequencies and the corre- 
sponding eigenvectors. The three frequencies are [11], 
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where in the first equation the plus and minus signs corre- 
spond to axial and radial modes with frequencies fl z = u>i 
and £l r = 0J2 respectively. The third frequency (which is 
the same as one of the frequencies in the case of the spher- 
ical trap) is the frequency of the two remaining degener- 
ate modes with Z = 2,m = ±l±2. For the cigar-shaped 
traps, i.e. in the limit of a very small A, the axial and 
radial frequencies reduce to 
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For the equation of state with the polytropic index 7 = 
2/3. the axial and radial frequencies are, respectively, 
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The corresponding velocities, Eq. (11), in the limit as 
A — > are given for the axial mode by, 
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and for the radial mode by, 
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where 7 = 2/3 was used. 

The frequency of the axial mode in Eq. (14), agrees 
very well with the value of approximately 1.55 measured 
by Bartenstein et al. [4] The frequency of the radial mode 
at the Feshbach resonance measured by Kinast et al. 
(fi r /uv = 1.829) also agrees remarkably well with the 
value given in Eq. (14) [3]. 

However, the frequency of the radial mode measured in 
the experiment by Bartenstein et al. ( [4]) is 1.67 which 
is about 20% below the corresponding value in Eq. (14). 
This deviation complements a sudden jump in the fre- 
quency and especially in the damping rate of the radial 
mode [4] not far from the resonance on the "BCS" side. 
There is no jump at this particular value of the scatter- 
ing length in the axial mode, so this phenomenon cannot 
possibly be associated with a change in the equation of 
state or transport properties by itself. Thus, an extra 
source of dissipation in the radial mode in the experi- 
ment by Bartenstein et al. must thus be associated with 
a non hydrodynamical effects 5 . 

We now apply viscous hydrodynamics to describe the 
damping of the collective modes in the vicinity of the 
Feshbach resonance. The primary source of dissipation 
in the hydrodynamic limit is shear viscous flow [12]. The 
rate of change of the total energy of a mode is given by 
[9] 



dE 
~dt 



= - 77 a 

3 



f r) ( dvi dv k 2 \ 2 3 

-J)' / r, ( r)Sr, ,17) 



where r\ is the coefficient of the shear viscosity and the 
ratio a z /a r is given in Eqs. (16) for each mode. 



^Such traps with A « 1 are referred to as cigar-shaped or 
prolate. 

^Such flow with velocity components proportional to the po- 
sition is often referred to as Hubble flow. 



5 We are grateful to R. Grimm for a discussion about a pos- 
sible source of this effect. 
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A damping rate of a collective mode can now be ob- 
tained by dividing the rate of change of energy, Eq. (17) 
by half of the the time-averaged total energy of the vibra- 
tional mode which is equal to maximum kinetic energy 
of the mode 
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where Eq. (3) was used with 7 = 2/3. In the last step the 
product of Thomas-Fermi radii was expressed in terms of 
the total number of particles N in the cloud, 



N = J n eq (r)d 3 r = J n eq (0) ( 1 " £ 
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and R x = R y = R r for the axially symmetric trap. 
Using Eqs. (17) and (18) the damping rate is 



(dE/dt) 



(E) 
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This expression will be used in the next section to extract 
the coefficient of shear viscosity of the strongly coupled 
atoms near the Feshbach resonance. 



III. QUANTUM VISCOSITY 

A. Generalities 

In the weak coupling regime of a small scattering 
length a, transport properties are determined by the par- 
ticle mean-free path 
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In the strong coupling regime, when a — > ±00, Eq. (21) 
becomes meaningless. Since the cross section does not 
diverge but is bounded by unitarity, this regime is some- 
times referred to as "unitarity limited" one. Indeed, 
the maximal possible two-body cross section is limited 
by Cmax = 47r/fc 2 for fixed collision energy (k is the 
wavenumber of relative motion). Thus, one may think 
that the mean free path is actually l m i n = l/w max . 
However, this is too naive as well, since at strong cou- 
pling there is no reason to limit kinetics to a picture of 



propagating particles rarely suffering only binary colli- 
sions. 

Whatever the microscopic structure of matter may be, 
the viscous hydrodynamics is an adequate description of 
low frequency dynamics. It is important that it is based 
on expansion in inverse powers of the cross sections, or 
expansion in small mean free paths I, so the stronger 
the interaction the better this approach is expected to 
work. The viscosity is in general defined via the dissipa- 
tive part of the stress tensor and can be defined without 
any assumption on the underlying matter. It appears in 
observables like the sound dispersion law 
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and thus can be measured. In fact, the damping rate in 
eq. (20) is the corresponding analog for trapped atoms. 



B. Universality of transport coefficients 

For simplicity we discuss mostly the zero temperature 
limit in this section. The temperatures in the experi- 
ments of of Kinast et al. [2,3] and Bartenstein et al. [4] 
range from 0.1 T F to 0.03 T F . 

One of the basics of low temperature physics of weakly 
coupled systems is that at low T the mean free path 
of quasiparticles goes to infinity, together with viscosity. 
There are multiple examples of such behavior, e.g. in 
liquid 4 He the minimal viscosity is around T c and then it 
rises as T — > 0. The normal component at low T is thus 
very viscous. 

However, this picture is no longer true at weak cou- 
pling. Strongly coupled normal component is also ex- 
pected to have low viscosity even at low T **. Further- 
more, we suggest that such oscillations even at T — > 
would have a damping described by a "quantum vis- 
cosity" proportional to the Plank constant, reflecting 
limitation on particle localization even at zero temper- 
ature. Indeed, at zero temperature a single interaction 
parameter — the scattering length a — diverges at the Fes- 
hbach resonance. The only remaining length scale is 
given by the inter-particle distance, n -1 / 3 . Thus, the 
mean free path must be of this scale. 

These arguments lead to universal relations for trans- 
port coefficients. In particular, the only form of the vis- 
cosity is 



77 = nha v 



(23) 



"For clarity, we remind again that experimentally observed 
oscillations are not an elementary quantum phonon state, 
but a highly excited state to which we apply the viscous 
hydrodynamics. 
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where an universal dimensionless coefficient a v is intro- 
duced. Note that it is independent on the particle mass. 
The same arguments are also valid for strongly coupled 
bosonic atoms, although of course with a different a v . 
At nonzero temperature viscosity has a similar form, 



77 = nha v + s(T)hj3 v 



(24) 



where the second temperature-dependent term contains 
the entropy density and another universal coefficient (3^. 
This second term reproduces the finite temperature, de- 
pendence of the damping observed by the Duke group, 
although we do not have any other justification for this 
form of this term. We repeat, that none of these terms 
are reproduced by binary collisions, neither in magnitude 
nor in parametric dependences. 

Similarly, the characteristic time scale in strongly cou- 
pled liquid is given by the Fermi energy and the Plank 
constant h, r _1 <~ ep/h <~ hn 2 / 3 /m. even if it in- 
clude purely Bosonic atoms. Thus, the scattering rate 
should be simply proportional to this scale with another 
universal coefficients, although different for bosons and 
fermions. 

The next question is what are the values of these di- 
mensionless parameters. For an ordinary liquid, these 
constants are of order one. For a quantum liquid, such as 
4 He with a large quasiparticlc mean free path the value 
of a v exceeds that of ordinary liquids by about three or- 
ders of magnitude. 

If the dimensionless parameters are much smaller than 
one, we deal with exotic near-ideal liquid. Although the 
available data available is not conclusive. It is is quite 
possible that the ultracold fermionic atoms near the Fes- 
hbach resonance are in the regime of near-ideal liquid. 



C. Experimental estimate for universal viscosity 
coefficient 

It was shown in section II that the damping rate of the 
small collective oscillations is proportional to the volume- 
integrated viscosity, Eq. (20). The universality relation 
Eq. (23) reduces the integral in Eq. (20) to the total 
number of particles times the coefficient a v which we 
want to determine. Thus, the damping rate is tt, 
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For the axial mode in the cigar-shaped potential trap 
using the ratio a z /a r from Eq. (15), the coefficient a v is 



a r , 



2hmR 2 z T z 
1927i 



(26) 



The experimental values are T z /uj z — 0.0036 at the point 
of the Feshbach resonance, with the minimal value of 
(T z /u z )\min = 0.0015 slightly off the resonance [4]. Using 
the axial trap frequency uo z ~ 140 Hz we get our final 
result for the minimal dimensionless viscosity- 1 - * 



a 



Ml 



0.3. 



(27) 



The fact that we get a number of order unity is an in- 
dication that the suggested "quantum viscosity" indeed 
exists and is described by the universality arguments. 

For the radial mode, using Eq. (16) and the fact that 
R z = R r /X we get, 



Otr, 



3m R 2 x T r 
8h 



(28) 



In the experiment of Bartenstein et al. [4] with the 
radial trap frequency uj r w 4700 Hz the damping rate 
of the radial mode at the Feshbach resonance is r r = 
0.0625 Lo r . Thus 
damping rate of the radial is 



the coefficient a v extracted from the 
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(r) 



1.1. 



(29) 



The coefficient a v extracted from the damping rate of 
the radial mode measured in the experiment of Kinast et 
al. [3] is 



0.2. 



(30) 



The damping rate of the radial mode measured by Ki- 
nast et al. [3] at temperature comparable to the Fermi 
energy is more or less consistent with an estimate from 
kinetic theory discussed in the next section. 

Kinast et al. [3] measured the temperature dependence 
of the damping rate. As the temperature is lowered the 
damping rate decreases in contradiction to the kinetic es- 
timates in section IV, but in agreement with our reason- 
ing. Unfortunately, the temperature is not small enough 
to conclude whether the universal "quantum viscosity" 
regime at T = is reached and a n 7^ or not, and 
whether the same universal viscosity is responsible for 
damping of both axial and radial modes. The second 
term in Eq. (24) proportional to the entropy density can 
explain this data if value of /3 V is of order unity. 



"The universal relation should not be valid near the edges 
of the system, at less than one mean-free-path or at optical 
depth less than one, where dissipation is larger. However, the 
edge includes only about one percent of particles and their 
contribution can be neglected. 



"The error is comparable to the value itself, as can be seen 
from experimental data points. Ironically, the situation with 
dimensionless viscosity of quark-gluon plasma is quite similar. 
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IV. COMPARISON WITH THE TRADITIONAL 
KINETIC THEORY 

In this section we will show that the kinetic theory 
based on the notion of binary collisions fails to describe 
the damping of the axial mode discussed above. 

Let us start with an order of magnitude estimates of 
the collision rates and viscosity neglecting Pauli block- 
ing and using the largest "unitary limited" cross section 
a = An/hp. The collision rate at the center of the trap 
estimated like this gives 



'coll 



n(0)avf 



10 5 £ 



(31) 



where the last number corresponds to conditions of the 
experiment of Bartenstein et al. [4] . Compared with os- 
cillation frequencies, w r = 4712 Hz, u) z = 141.9-ffz of the 
trap, which leads to a conclusion that only the latter 
mode have a chance to be hydrodynamical. 

The mean free path, l m i p , of a particle is of order 
l/(ner), while the shear viscosity is 



_ , mv 

r\ ~ mvn W P = , 
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(32) 



where v is the average velocity of a particle. In the 
limit of zero temperature the velocity is set by the Fermi 
momentum, mvF — hkp. In the vicinity of the Fes- 
hbach resonance the cross section is unitary bounded, 
a < o max = An/kp. So, if we take its maximal value 
(and still ignore for a moment Pauli blocking), we will 
get a minimal viscosity which may follow from binary 
collisions: 



n 40 
hn 6tt 



(33) 



This inequality is strongly violated in experiment, as 
shown above: this "minimal binary" value is in fact larger 
than observed value, and forty times larger than the de- 
rived bound possible for a liquid. 

Furthermore, if the fermionic atoms were in the Fermi 
liquid regime, the collision rate will be significantly low- 
ered by Pauli blocking which should lead to a suppression 
factor about (T/Tf) ~ 1/1000 in the experimental con- 
ditions. If true, the oscillations then would be basically 
collisionless and no hydro phenomena would be present. 
In a picture of BCS-type pairing, with relatively small 
modification of Fermi sphere, T in the above formula 
to be substituted by a gap, so the rescattering suppres- 
sion would be of the order of (A/ep) 2 - The gap value is 
not well known, but this suppression factor is still about 
1/100 or so. We must then conclude that both pictures 
are wrong and in fact there seem to be no Pauli blocking 
whatsoever^. 



To make more quantitative conclusion we will derive 
here the damping rate of a collective mode in a axially 
symmetric trap applying the traditional kinetic equation 
to an almost ideal Fermi gas with unitary limited cross 
section. 

A damping rate of a collective mode in the kinetic the- 
ory is determined by a relaxation time which is a mea- 
sure of how fast a particle distribution function n(p, f, t) 
for a given collective mode takes an equilibrium form. 
Both the time dependent and equilibrium distribution 
functions are the solutions of the kinetic equation. The 
equilibrium distribution for a Fermi gas is 



e (e- M +V(r))/fc B T + 1 



(34) 



where e = p 2 /2m. 

During an oscillation the distribution function is dif- 
ferent from the equilibrium one. The collisions between 
particles cause the non-equilibrium distribution function 
to "relax" to the equilibrium form. These collisions are 
the source of the damping of the oscillations. 

As shown in [10] the damping rate of the oscillations 
of Fermi gases is equal to 
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where 
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(36) 



(37) 



where A$ = ($i + $ 2 ) - (<& 2 + $2) an d the function 
<f> = p 2 — ^j- describes the deviation of the distribution 
function of a collective mode from the equilibrium one. 
The function W is proportional to the scattering ampli- 
tude of binary collisions pi , p 2 — > p^ x , j? 2 ■ ^ n tn c vicinity 
of the Feshbach resonance this function is determined by 
unitary limit of the scattering amplitude and is equal to 



§§ After observation of elliptic flow this issue was discussed 



in literature and the MIT group [13] have argued that this 
might have been due to strong deformation of Fermi sphere 
in exploding gas. This explanation obviously would not be 
applicable to small amplitude oscillation we study here. 
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p = pi — P2 is the relative momenta of two particles. 

At a very low temperatures the Pauli blocking factors, 
n (1 — n), in Eqs. (36) and (37) significantly reduce the 
phase space of particles whose collisions appreciably con- 
tribute to the relaxation. The main contribution is from 
the collisions of particles whose momenta lies very close 
to the Fermi surface: 



e-fji + V(r) ~k B T 



(39) 



where fi is the chemical potential Eq. (9). 

After a lengthy but a straightforward calculation one 
obtains for a damping rate: 



97T (k B TY 
~ 50 hfi 



(40) 



It is important to note a temperature dependence of the 
damping rate. It has a typical T 2 dependence which 
comes from the life time of weekly interacting Fermi par- 
ticles or quasi-particlcs in the case of the Fermi liquid 
[9]. We stress that such scaling is true only for weekly 
interacting gas of particles or quasi particles which is ob- 
served in liquid 3 He. If the atomic cloud near Feshbach 
resonance is indeed a strongly interacting near perfect 
liquid the damping rate and other dissipative processes 
will very weekly depend on the temperature. This pre- 
diction can be checked experimentally. 

The damping rate can be expressed in terms of the 
trap frequencies and the Fermi temperature Tp as 



r,| (3W ) 



1/3 _ / t y 



(41) 



For the experiment of Bartenstein et al. [4] the above 
formula predicts for the axial mode 



— w 0.56. 

ui z 



(42) 



This is much larger than the damping observed close to 
the Feshbach resonance, and is only compatible with data 
well away from it, where Uf a < 1. 

Similarly, for the experiment of Kinast et al. [3] the 
Eq. (41) predicts the ratio for the damping of the radial 
mode to be 



— w 0.17, 



(43) 



which is almost ten times larger then the observed ratio 
T r /uj r w 0.014. 

Thus, the strongly interacting atomic cloud is better 
described by a picture of a nearly perfect strongly in- 
teracting liquid with viscosity very close to the minimum 
possible value. As one moves away from the Feshbach res- 
onance into regime of weekly interacting Fermi gas the 



kinetic theory becomes again applicable. Furthermore, 
as detuning gets larger than used in the experiment dis- 
cussed, the scattering length is getting small, the gas 
will enter an almost collisionless regime, and the damp- 
ing rates is getting small again. Thus the damping rate is 
expected to reach a maximum value for certain value of 
the magnetic field, separating weakly and strongly cou- 
pled regimes. 

For the radial mode the value predicted by the same 
kinetic calculation is 



— w 0.02. 

UJ r 



(44) 



comparable to the value of about 0.06 obtained experi- 
mentally. It shows that this mode may not be treated by 
hydrodynamics. 



V. DISCUSSION 
A. Other strongly coupled systems 

In the introduction we mentioned a number of strongly 
coupled systems which display similar transport proper- 
ties. One example is the strongly coupled quark-gluon 
plasma (sQGP) [5] which was found in heavy ion exper- 
iments at RHIC at temperatures above the critical tem- 
perature, T = (1 — 2)T C = 170 - 350 MeV. In ultrarela- 
tivistic heavy ion collisions also quite spectacular explo- 
sions are observed, with radial and elliptic flows surpris- 
ingly well described by ideal hydrodynamics. The sQGP 
seems to have a very small viscosity, i.e. 77/5 « 0.1 — 0.2 
[14,15] and rj/s w 0.2 — 0.4 from lattice simulations [16]. 
It is not even far from the lowest limit discussed below. 
Two of us suggested that the sQGP has such low viscos- 
ity because of the existence of weakly bound states near 
the so called "zero binding lines" on the QCD phase dia- 
gram [5] . It was gratifying to learn after that that similar 
role for trapped atoms is played by Feshbach resonances. 

Other examples are: strongly coupled J\f = 4 Super- 
symmetric Yang-Mills (SYM) gauge theory, a four di- 
mensional conformal field theory (CFT) and the strongly 
coupled QED plasma. The classical one-component QED 
plasma is particularly well studied and is known at val- 
ues of the effective coupling T = (Zefn 1 ^ /T = 2-300 
to be a liquid, with a very small viscosity which has its 
minimum at T ~ 10 (see e.g. [17]). 

The CFT is a toy model, emerged mostly in the con- 
text of string theories. It is a model resembling QCD, 
the gauge theory of strong interaction. In CFT as op- 
posed to in QCD the gauge coupling is not "running" 
with energy scale. Thus, the strong coupling regime is 
obtained simply by using a large coupling constant in the 
Lagrangian, A = g 2 N. (Here g is the gauge coupling and 
N is the number of colors. Only this combination ap- 
pears if N is large.) This limit can be addressed using 
the AdS/CFT correspondence as originally suggested by 
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Maldacena [18], whereby the quantum intricacies of the 
strongly coupled gauge theory are mapped onto a classi- 
cal problem in gravity albeit in ten dimensions. The finite 
temperature version of this theory describes a plasma- 
like phase with strongly coupled constituents. The four 
dimensional world (in which the CFT fields live) is a sur- 
face in ten dimensional space, at some distance from a 
black hole, with a mass adjusted to yield the desired tem- 
perature T at this surface. We would like to mention two 
important results that follow from this construction. One 
is the equation of state of the underlying gauge theory at 
strong coupling A> 1 [19] 



PA CO 

Po(T) 



1 



= (l-^ + 0(l/A 3 / 2 ) 



(45) 



where po (T) ~ T 4 is the Stephan-Boltzmann pressure for 
zero coupling. The second result is the viscosity of the 
underlying gauge theory at strong coupling [20] 



X^ = i( 1 + °< 1 ^ /2 >) 



(46) 



given in units of the free entropy density s***. The cor- 
rections were recently calculated in [21]. 

While the pressure is only changed by 1/4 when the 
coupling is changed from zero to infinity, the viscosity 77/s 
changes from infinity to a finite (and surprisingly small) 
number. Thus, one may wonder whether other strongly 
coupled systems show similar behavior, and whether such 
limiting numbers can be universal and theoretically un- 
derstood. 

The holographic principle in the Maldacena limit and 
the Kubo formulae show that the viscosity is proportional 
to the graviton absorption cross section in bulk by the 
black hole, while (according to Bcckcinstcin-Hawkins ar- 
gument) the free entropy is related to its area. As a re- 
sult the same limit for the viscosity holds for a number of 
backgrounds, even in different space dimensions. These 
observations led to a conjecture that Eq. (46) is a univer- 
sal lower bound valid for any thermal system in strong 
coupling [22] . Below we provide its heuristic derivation of 
based on the uncertainty relation and Einstein's famous 
relation between the diffusion constant from the fluid vis- 
cosity. As a result, we show how the bound (46) fits well 
into the the liquid-like picture of CFT at finite tempera- 
ture. We then use these insights to derive an even lower 
bound for cold Fermi systems, and conjecture that at 
strong coupling they also make an universal near-ideal 
liquids. 



B. Bounds on Transport Coefficients 

For atomic systems one may also think that as the 
interaction strength is driven to infinity the transport 
parameters such as viscosity become as small as possible. 

In a weak coupling limit the viscosity and diffusion co- 
efficients are both related to the scattering length and are 
thought of as proportional to each other. For liquids one 
should think differently, An example of quite an opposite 
relation between them was provided by the famous Ein- 
stein relation, which we now derive for consistency of the 
presentation. 

The distribution of suspended particles in a thermal- 
ized column of gas is given by statistical mechanics. In- 
deed, if n(x) is the suspension density at finite tempera- 
ture then 



n(x) 
MO) 



-mgx/ksT 



(47) 



which follows from Boltzmann. Einstein's observed that 
an arbitrary sphere of radius ro in suspension within the 
column would also follows the same "distribution" pro- 
file. The idea is that the sphere under gravity will fall 
with a terminal Stokes velocity 



vt 



mg 
6ir r rj 



(48) 



but the fall will be balanced by random upward kicks 
due to Brownian motion. In equilibrium, the upward 
diffusion balances the downward gravitational fall so that 
in the stationary limit 



dx 



-hvt 



(49) 



from which it follows that n(x)/n(0) = e ~ v Tx/D^ q 0U1 _ 
paring this result with Eq. (47) and using Eq. (48) yields 
the Einstein's formulae 



D 



k B T 
6ir ro T] 



(50) 



Although this formulae was derived for a macroscopic 
sphere of radius r immersed in a suspension, empiri- 
cally it is known to hold through 14 orders of magnitude 
changes down to the suspension constituent wavelength 
[23]. 

We recall that in three dimensions the diffusion con- 
stant is just D = v 2 t/3 = A 2 /(3t) where A and r are 
the mean- free path and collision time ttt. Inserting this 
result into Eq. (50) yields 



**Tn thermal gauge theories as in the case of blackbody 
radiation, there is no ordinary particle density but only the 
entropy density s ~ T 3 . 



tt1 We note that in p space-dimensions the diffusion con- 
stant is D = A 2 jpr and (50) should be derived accordingly. 
All the bounds to follow can be extended readily to p space 
dimensions. 
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v(feri (fcs7V) 



(51) 



In a densely packed liquid the smallest jump (the 
mean- free path A) is the size of the quasiparticles r . 
Classically in densely packed hard balls t can be as small 
as zero due to the fact that they are always touching. 
Quantum mechanically however this is not allowed since 
the time localization cannot be better than the limit set 
by the largest allowed energy, by the Heisenberg uncer- 
tainty principle, i.e. ksTr > Ti/2. Inserting this result 
into Eq. (51) yields Eq. (46) since the entropy per unit 
volume s is just the number of (quasi)particles per unit 
volume due to the close packing, i.e. s/ks = n = 1/A 3 . 
The ensuing physical picture of the strongly coupled ther- 
mal system is that of a liquid with the shortest time cor- 
relation length T min = h/2kB T. 

Our heuristic derivation of Eq. (46) follow from the 
assumption that Eq. (50) holds for the liquid particles, 
since the relation is known to hold over many orders of 
magnitude changes in n, D. Thus the particle and en- 
tropy densities are the same. While classically the colli- 
sion time is zero for the densely packed liquid, quantum 
mechanically it is bounded from below by the Heisenberg 
uncertainty principle. Thus, 



(52) 



which is the same as the CFT limiting value. Turning 
the argument around through Eq. (50) implies an upper 
bound on the diffusion constant in a strongly coupled 
liquid, namely 



D < k B T 

cr ~ 12 ft 



(53) 



in three dimensions, where the cross section cr = 87rrg. 

Let us now turn to cold atomic gases and repeat the 
same argument once more. It is simpler to imagine a 
Fermi gas in a vertical gravity field, for which we will re- 
run Einstein's derivation. We note that the trap field ac- 
tually fulfills the same role. Using either Thomas-Fermi 
or hydrostatic calculations one finds that the Fermi mo- 
mentum for non-relativistic quasiparticles in a gravita- 
tional field is 



ft&F = (2m* (fj, — mgx)) 



1/2 



(54) 



where m* is the quasiparticle mass and m is its bare 
gravitational mass. The normal Fermi density in a weak 
gravitational field is 



n(x) f mgx\ 



n(0) 



3/2 



Smgx J (2fx) 



i 1 ) 



(55) 



Using Eq. (55) and repeating steps which led to Eq. (51) 
one gets, 



n 



V(r A 2 ) 3tt 



(56) 



with a diffusion constant 



D 



9tt ro?? 



(57) 



In the infinite coupling limit A becomes ro and the sys- 
tem is again closely packed. The Heisenberg uncer- 
tainty principle stipulates that the shortest collision time 
is dictated by the largest available quasiparticle energy, 
namely [mt > ft/2. Thus the new bound on the viscosity 



Hn n ~ Qtt 



(58) 



Although Eq. (58) was derived for non-relativistic par- 
ticles, its insensitivity to the quasiparticle velocities im- 
ply that it should hold in the relativistic case as well. 
Equation (58) implies an upper bound on the diffusion 
constant 

in trapped cold Fermions in three dimensions. 



VI. CONCLUSIONS 

Strongly interacting systems, in relativistic field the- 
ories (QCD, CFL) and in condensed matter physics 
(strongly coupled QED plasma, ultracold atoms near the 
Feshbach resonance), are radically different from gas- like 
weakly interacting systems. The best way to see that is 
not via the equation of state and related thermodynami- 
cal observables, but with the help of transport properties. 

As we have argued above, all such systems are near- 
ideal liquids, and thus the natural tool one should use 
to describe those are standard viscous hydrodynamics. 
Indeed, one gets very good description of the "elliptic 
flow" , the frequency and damping of the lowest collective 
oscillations. 

Furthermore, the viscosity extracted from the data are 
shown to be very different from what is expected on 
the basis of binary collisions and in the weak coupling 
regime. We compared the minimal experimental value 
of the viscosity (or maximal rescattering rate) Eq. (27), 
with results of standard kinetic theory. Even without 
Pauli blocking and with maximal (unitary limited) bi- 
nary cross section the kinetic theory fails to reproduce 
the data. We have also shown that the kinetic theory de- 
scribes damping far from Feshbach resonance, where the 
the weakly coupled gas-like regime is valid. We conclude 
that the system of optically trapped ultracold atoms near 
the Feshbach resonance, like other strongly coupled sys- 
tems discussed in section V, is not a gas but rather a 
near-ideal liquid. Since oscillations are maximal at the 
surface, this small viscosity does not refer to superfluid- 
ity and is property of the normal component of the liquid 
making. Thus, the the properties of the normal compo- 
nent of the atomic system near the Feschbach resonance 
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is quite different from that of liquid 4 iJe below the A- 
point. 

At infinite coupling the constituents are effectively 
large and densely packed. The packing fluctuates over 
short time scales dictated by the Heisenberg uncertainty 
principle. These physical insights together with Ein- 
stein's description of diffusion in viscous liquids, allow 
for a simple derivation of the viscosity to entropy den- 
sity ratio established using CFT. We have extended this 
derivation to cold Fermi systems and derived an even 
lower bound for the viscosity to the particle density ra- 
tio. 

It would be quite important to make more accurate 
measurements of the damping rate in order to see how 
close is the minimal viscosity to the theoretical bound. 
The temperature dependence of the damping rate would 
clarify the issue of Pauli blocking. We expect that the 
temperature dependence should be weak and not even 
close to what binary collision theory for Fermi gas pre- 
dicts: the system is not even qualitatively close to a Fermi 
gas at the strong coupling limit. 

Temperature dependence is also crucial for under- 
standing of the transition to superfluidity which was 
avoided in this paper. We argued that in super- 
component resigns in the interior of the system, which is 
not very important for elliptic flow or oscillations we dis- 
cussed. It seems that the lowest oscillation mode is well 
described by the one-component hydrodynamics: how- 
ever it is of course quite likely that some higher excita- 
tions are analogous not to the usual sound but to other 
sounds known for supcrfluid liquid He. 

One should also develop a theory of the strongly cou- 
pled systems beyond quite schematic mixtures of the BCS 
superconductor and ideal Bose gas of molecules. So far, 
such models can approximately reproduce the equation 
of state but not the transport properties. Numerical sim- 
ulations of larger scales can also be helpful: perhaps one 
should complement the equation of state calculations by 
measurements of long-time correlators related by Kubo 
relations to transport coefficients. 
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